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TECHNICAL NOTE NO. I833 


NOTES ON THE FOUNDATIONS OF THE THEORY OF SMALL 
DISPLACEMENTS CF ORTHOTROPIC SHELLS 
By F. B. Hildebrand, E. Reissner, and G-. B. Thomas 

SUMMARY 


A survey has "been made of various systems of equations ■which, have 
been given in the literature for' the analysis of small deflections of 
thin elastic shells. In this survey the results, previously known for 
isotropic shells, have been reformulated for shells which are orthotropic 
to the extent that the normal to the middle surface of the shell may be 
a preferred elastic axis. 

A new system of equations has been derived for the analysis of 
shells, which includes the effects of transverse shear and normal 
stresses. The assumed orthotropy of the shell facilitates the 
identification of the separate effects of the ordinarily neglected 
transverse stresses. 


INTRODUCTION 


In this report various methods are described of obtaining systems 
of equations which may be considered as forming the basis of a theory 
of small displacements of elastic plates and shells. In a.1 1 cases 
one begins with the governing equations in the three-dimensional 
theory of elasticity and an attempt is made to reduce this system of 
equations. Involving three Independent space variables, to a new 
system involving only two space variables. These two variables are 
most conveniently taken as coordinates on the middle surface of the 
plate or shell. 

Numerous reductions of this sort have been carried out in the past 
by different workers in the field of elasticity. Care will be taken to 
point out in the analysis the connection between the developments herein 
presented and the results of earlier studies. 

One of the points of interest is the study of the effect of 
transverse shear deformation on the bending of Bhells. In the case 
of the flat plate It has been Shown, by methods which differ somewhat 
from those used here (references 1 and 2), that inclusion of this 
effect resolves in a natural way well-known difficulties with regard 
to the boundary conditions which may be prescribed along the edges of 
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a plate. • Itr^La to “be expected that corresponding difficulties in the 
analysis of shells may be removed in an analogous way, and. the present 
report— shows that tills is indeed- the case. The practical significance 
of- this analysis with reference to stress— concentration problems has 
been established in the earlier work on plates. 

Attention is herein restricted to static problems involving small 
deformations, excluding for the present the study of vibrations,' elastic 
stability, and finite deformations. It will be apparent, however, that 
much of the present analysis can be extended to the consideration of 
these problems. 


In order that the separate effects of the transverse stresses may 
be segregated, the plates and shells are assumed to consist of material 
which may be orthotropic to the extent that the normal, to the middle 
surface is a preferred elastic axis. 

This work waB conducted at the Massachusetts Institute of Technology 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 


Si* S2 

e 

A l* a l* “2 

E l* ^2 
XT 

! 

€ 1 * e 2> 

7 12> 7 l£* 7 2 £ 
ff l, a 2* 


coordinates in middle surface of shell 

coordinate normal to middle surface 

unit vectors in directions of £g, £ 

parameters in linear element (equations (2) 
and (3)) 

principal radii of curvature 
displacement vector (jj - Uj_tj_ + TJgtg 
body— force vector ^F = F-jt-^ + Sgtg + 
components of direct strain 
components of shearing strain 

components of direct stress 


+ Wn) 

y c H ) 
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t 12j t 1£> t 2£ 
V, V g 

G, 

V 

V* 

E 

it 

P 

P 

m 


n l , n 2 


Mj2 


P ll* P 12-» P 21-» 

Si, S 2j T]_, 

*2 

A, B 


u l> u iS u i M 1 

u 2 , u^*, u 2 ” r 

V, w’ J w” 

J 

0 J tt 1 

i 

? , e , 6 


3°j P\ P" ! 

>* 

n°, n* s u"j 



components . of shearing stress 

Poisson’s ratios 
Young’ s moduli 
shear moduli 

parameter defined by v = 


E 


parameter defined "by v* - 
parameter defined "by E = - 


VfE 


(1 - v)E^ 
1 - V 


E 


) 


1 _ v - 2v 

potential energy 

strain energy per unit volume 

effective external force per unit area; 

applied to middle surface (p =* p-jt-^ 4 - Pgtg + qn) 

effective external moment per unit area; 

applied to middle surface (m = — mitpr) 

streas— resultant vectors 

(^1 = ^ll^l + ^12^2 + ^l s * % = ^21^1 + N 22 t 2 + 

stress— couple vectors 

( M 1 “ M 12 t l ~ M ll t 2 5 ^2 = M 22 t l ” M 21 t 2) 

higher-order stress resultants defined hy 
equations (73) 

auxiliary stress resultants defined "by 
equations (79) 


displacement functions defined hy 
equations ( 69 ) 


strain functions defined hy equations (8 ^) } 
■with appropriate subscripts 
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7°, k, t strain functions defined "by equations ( 46 ) 

and ( 48 ft), with appropriate subscripts 

a 0 , a*, a ft auxiliary parameters defined "by equations (86), 

with appropriate subscripts 

THE?EB- 3 )IMEKSOTAL STRESS AMD STRAIN' IN CUHVTLINEAR COORDINATES 


In this section there are collected for convenient reference 
certain taiown ftasic formulas pertaining to the analysis, of stress and 
strain In terms of orthogonal curvilinear coordinates. 

THE COORDINATE SYSTEM 


The position of any point in a plate or shell may he specified by 
three coordinates S 2 , £ * where and £2 s-peoify position on 

the middle surface, while £ measures the distance, along the outward 
normal, from the middle surface to the point. In order that the 
coordinate curves be orthogonal. It is required that the £•]— 

and Sg— curves be the lines of curvature on the middle surface £ = O. 

The unit normal vector at a point of the middle surface is denoted 
by n and the unit tangent vectors to the and ^2“ cur ves are 

denoted by and Tg, respectively. The coordinates and | 2 

are to be chosen in such a way that the system is right-handed , in the 
sense that ¥1 Is rotated into t 2 by a right-handed rotation 
about n. 

If there Is written for the position vector to "appoint in space 

E (£i* hj £) • *2) + £ 5 ( s i^ £2) Ci) 

where r Is the position vector to a point on the middle surface, the 
linear element Is of the form 


ds 2 = cLR*cLR = A-] 2 d 5 x 2 + A 2 2 d£ 2 2 + d£ 2 


( 2 ) 


The coefficients in the linear element are given by 
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where H-j_ and Kj at*© the principal radii of curvature of the middle 
Burface and 


_ 2 _ Sr Sr 
1 ^1 a &l 


1 _ 1 Sn _ Sr 

aj^ S^ 


„ 2 _ Sr Sr 
2 Si 2 S| 2 

1 _ 1 Sn i Sr 
R 2 a ^ 2 Sl 2 S ^ 2 


(M 


There is also the relation 


n = x t 2 


L_ Sr. . * Sr 

^2 agi ss 2 


(5) 


The further relations 



(6) 


which are special cases of formulas obtained by Mainardi and Codazzi, 
are of frequent use. It is noted^ for later reference s that 
equations ( 6 ) imply the truth' of the equations 



( 7 ) 


COMPONENTS OF STRAIN 


If the displacement vector IJ is written in the form 
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the ala: components of -Infinitesimal strain are expressed in terms of 
the components of Tf by the equations 

(T^ 



y = 5 r u g 1 + + sAi) 5 r pj 

12 “l^ + ^1 a 2(-*- + £^2) + £^2) ^2 + £/%) 


+ [1 + 


i0^(l7 


’’as - 


02(1 + £/%) ^2 


+ ^2 


DIFFERENTIAL EQUATIONS OF EQUILIBRIUM 
If the body force per unit Yolume is denoted by the rector 


F = F-jt^ + Fgtg + F^n (10) 

the six components of stress 03 ., 02 , a £y r 12y t 1^j an ^ T 2£ must satisfy 
the three equilibrium equations 


r f 



NACA TN No. I833 


7 


f. 


\ 


^-[og(l + S/B^oj.] + + 5 /R!)t 12 ] 

+ a l a 2 [(i + + f/RsjT^] 

+ Tl2 + S/E i)] + T i£ ^0- + 5/%) 

“ a s Se^faC 1 + £^s)] + °i“2(i + £/Ri)(i + e/Rs)*! = 0 


^K 1 + ^ 1 )° a ] + af^C 1 + ? ^a) T ls] 

+ “aw fg-fC 1 + £Ai)(i + £Aa)Ta£] 

+ Tia si^hC 1 + £/Ba )] + T£ £ ^lirC 1 + ?/R i) 

i 

~ ai af^E^C 1 + ^1)] + + s/Ri)(i + ^2)^2 


Mil) 


0 


* 


a i°a 


|^[(i + S/Ri)(i + S/Ra)o s J + ^-[^2(1 + S/R 2 ) t i^ 

“ + 5/Ri)cr 2 + a ia g(l+ S/R^l + S/R 2 )^ = 0 



8 


NACA TN Ho. 1833 


STRESS— STRAIN" RELATIONS 


If It Is required that" the normal to the middle surface "be an 
axis of elastic symmetry and if isotropy is assumed in elements of 
the middle surface, the stress— strain equations may he written in 
the form 



*12 

7 lt 


. h 2 

G 

M 

* 

.Z& 




(12) 


(13) 


If equations ( 12 .) are solved fey and cfg the results can he 

expressed in the form 


E / \ V £ E 

ai ■ rr ^ 1 + V€2 ) + a- 

ji V f E 

cr 2 - U + veil + 7 — 

1 - v 2 ' V ' (1 " V ) E £ 



(HO 
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where 


a i = 


E£ 


(1 - V)E£ 

If there are Introduced the symbols 


V^E 


s + (1 - v)^ ( S 1 + s a) 


(14a) 


v * = 


't E 


E* = 


1 ~ V 

(1 - \> 2 )E£ 
1 - 2V^V* 


(15) 


equations (lh) and (lha) may also he written in the form 


°i - 7T^2 ( 6 i + ve a) + v *®£ 


°2 - 2 + V6 l) + ^ t 


e* r *. / \ 

^ p + v ( e i + 6 s) 


(16) 


An equivalent form of equations (lh) and (lha), which is more 
convenient for some purposes 3 is obtained in terms of the parameters 



( 17 ) 
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With tills notation aquations (l4) and (lta) lead to explicit 
expressions in terms of the Btrains in the form 


?! = ■ - ft ^2 (l - V )e x + (v + V )e 2 + V ft (l + v)<^J 
(l - V )e 2 + (v +■ V 2 )^ + (1 +• v)g£ j 




= E 


_rrr( s i + e a) + e {] 


08 ) 


(19) 


It is noted that in the limiting orthotropic case for vhich 


V* » 0 


(20) 


there follows V = 0, E = E, and equations (l8) and (19) "become 


1 - rr^ + v «e) 

2 - “^('2 + Ve i) 


a t = Es s 


( 21 ) 


Further, in the isotropic case for which 


V c = V 


E^ = E 


( 22 ) 
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\ 


there follows 


E 


1 - V 


isotropic ” d + v )d-2v) 


E 


v Isotropic 


= V 




(23) 


and equations (l8) and (19) "become 


- (1 + vK1-2V) [ (1 ~ V) ^ + + S 0] 


a 2 = 


E 


(1 + V) (1 - 2 V) 


|^(l — + 6 0j 


- V)e £ + V <U + '«$ 


(2k) 


THE MUTIMDM PRINCIPLE EQR THE DISPLACEMENTS; PRINCIPLE OE 
MENTMDM POTENTIAL ENERGY 


If no "body forces are acting and if all "boundary conditions are 
stress conditions, the principle sets forth that among all possible 
states of strain the state which actually exists is that one for which 
the potential energy ?r tales on its minimum value. (See, for example, 
reference 3j page 281.) The potential energy of the load system is 

given "by the surface integral 


*Z 



•U dS 


( 25 a) 


where p g is the surface-stress vector, while the energy Strain is 
given "by the volume integral 

^strain e £ J 7 12> 7 l£* 7 2fj) dT 
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where P is the strain energy per unit -volume which is to "be expressed 
in terms of the displacements Uq, and. V by use of equations (9). 

The function P is determined by the relations 


* 





dP 

^12 


T 12 


dP 


T i z 



(26) 


and, for- the orthotropic material considered here, takes the form 


p -§<! 


+ 


^~~ 2( € ! 2 + e 2 2 + ^*1*2) + G ^ 12 2 

+ v*( tl + < a )] 2 + o £ (r l£ 2 + r 2£ 2 ) 


(27a) 


Equation ( 27 a) can also be written in the form 




E 


1 — L 


(l - V 2 ) (.* + '*) + $(l - v2) e£ 2 + 2 (v + V 2 ) , x , 


+- 2 v £ (1 + v) ( e l + € 2) e £ + G 7 12 2 + ( 7 1 ^ 2 *~ 7 SZ/ 


(27b) 


The requirement that P be a positive definite quantity leads to 
the following restriction on the elastic constants 


E^ (1 - v) _ 2Ev^ 2 > 0 


(28a) 


/ 
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and hence the restrictions 


( 28 b) 


on the parameters V and V* or, equivalently, to the requirement 
that 1 and E* he positive. With the definitions of equations ( 25 a) 
and (25h), the minimal condition is of the form 

Stc = 5^ + Strain) = 0 (29) 


THE MIIOMJM ERINCIPIE FOR THE STRESSES] CAST IG-LIAHO 1 S THEOREM 

OF IEAST WORE 


If linear stress-strain relations and small strains are assumed, in 
accordance with equations (9 ) , and if the surface stresses are prescribed 
over the faces of the shell or plate, while the displacements are 
prescribed over the edge surfaces, the principle sets forth (see, for 
example, reference 3, page 286) that among all statically correct states 
of stress the true state is determined by the condition that the 
complementary energy jr c be a minimum 

8* c =0 (30) 

The complementary energy consists of two parts: (l) The volume 

integral of the strain energy, expressed in terms of the stresses, 
and (2) the work of the boundary stresses over that part of the 
surface over which the displacements are prescribed, 

* c = 1 ’ a 2 > t 12j T l£* t 2^) dV - JJ P-U 4S (31) 

edge 
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The function P Is determined, by the relations 


SP - 

ScTt ^ 


&p _ ~ 
^ T 12 12 


dP ' 

8 ffo 2 


_££_ = 
& T i5 


?ii 


gp-C 


(32) 


and, for the material under consideration here, is of the form 


r 


p = 

2 


Ci 2 + cr 2 2 + ff £ 2 ~ 2^1 °2 — 2 (cr-]_ + cr 2 )cr£ 


E £ 


+ 4 T 12 a + Oj( T 15 2 + T 2i S ) 


( 32 a) 


STRESS RESULTANTS AND COUPLES 


Equations ( 9 ), (ll), (13), (l 8 ),-~.and ( 19 ) comprise 15 equations 
involving the 15 unknown quantities Uj_, U 2 , W; ex, € 2 , , 7 x 2 j ?x£ j 

a l> a 2* r 12> T l£, t 2 $* quantities here depend upon the three 

independent-space variables |^, i 2 , and £ . 

For the purpose of obtaining a two-dimensional theory of plates 
and shells it is customary to eliminate the ^—coordinate in the 
expressions involving stresses by the Introduction of stress resultants 
and couples. Ten such quantities are conventionally defined, by the 
following equations: 
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N~1 1 = + £/& 2) ^ N 22 = I ^ 2(1 + ^/^l) ^ 

lJ-h/2 J-h/2 


N12 = 


/ t 

T-j^^ + tfaz) ^ ^21 = T l 2 (/- + £/%) ^ 

-h/2 (J-h/2 


ph/2 ph/2 

Ql = I T X^(^ + ^/^2^ 4 £ Q2 = T 2^(^- + ^ /^l} ^ > ( 33 ) 

ll-h/2 U-h/2 


lh/2 ph/2 

*11 = ^l(l + £/R 2 )£ ^ ^22 = ^(l + £/R;l)£ d£ 

U-h /2 U-h /2 


■\h/2 ph/2 

Ml2 - t 12 (i + £/R 2 )$ as %l - t i 2 ( 1 + £/Ri)5 ^ 

J-4i/2 U-h/2 


The significance of the 10 resultants and couples so defined is 
suggested by the laws of mechanics 3 irrespective of the material of the 
shell and irrespective of the state. of deformation of the elements of 
the shell. These resultants and couples can be considered as the 
components of the vectors 


N-j^ . = N-) 1 t-j + N-^gtg + Qjti Ng s= Ngitg 4- N pg tp 4- Qgn 


M 1 = %2^1 ~ ^2 


Mg = Mggt;^ - ^lt 2 


Tt is evident that these vectors represent resultant force and 
moment s per unit of length along the parametric curves 3 acting on 
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sections 1 -^ = Constant and = Constant. The couples are positive 

when they produce positive stresses on the part of the shell or plate 
on the positive side (£ > 0) of the middle surface. The absence of a 
third component in and is due to the fact that a differential 

element of area in a cross section has finite height h and infini- 
tesimal width. 

To obtain differential equations of equilibrium for the resultants, 
the equilibrium equations (ll) may be averaged over the thickness h of 
the plate or shell. Making use of equations (7) where necessary, the 
results are obtained in the form 


hshi + + Hl2 &2L _ Noo *2L + (h as + 


a Si 




8 I 2 22 81 ^ 


Qi 




a^Pi 


= 0 




(35) 


where 


Pi = 


+ JLYl +-i-Y 


Ei 


*2 ^ 


h/a 


H-h/2 


)h/2 


i-h/2 


1 + 4-Yi + d^ 


Ri 


P2 - 


1 + ifeX 1 


h/2 , 


q = 




’J-h/2 

h/2 


U 


-h/2 


ph/2 

-h/2 
h/2 

I-h/2 


1 + ^A 1 + feK 45 ^ (36) 


i + £ 


Q( x + 


a£ 
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The quantities p-^, Pg, and q can "be considered as the components of 
« - , a rector 

P = Pi*! + P 2 ^2 + ^ (37) 


representing effective external force per unit area applied to the middle 
surface of the plate or shell. 

Two differential equations of equilibrium for the couples are 
obtained by multiplying both sides of the first two of equations (ll) 
by £ and integrating over the thickness. It may be noted that if the 
third equation were treated in a similar way new quantities not defined 
by equations (33) would be introduced. The two equations described are 
obtained in the form 



The quantities and mg can be considered as the components of a 

vector 

m = mgt]_ — m^tg (40) 

representing effective external moment per unit area applied to the 
middle surface of the plate or shell. 
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Equations (35) and (38) comprise 5 equations involving the 
10 resultants and couples defined "by equations (33). Two different 
points of view may now "be adopted. 

In the first place, the 10 quantities defined by equations ( 33 ) 
may he considered as a complete and adequate macroscopic description 
of the statics of the plate- or shell. In this event, it becomes 
necessary to- determine which quantities serve to pr'ovide a 
corresponding macroscopic description of the state of deformation 
of the plate or shell which is logically equivalent to the foregoing 
description of the state of stress. By this is meant that in the 
description of the state of deformation there should appear no 
quantities whose values are affected by stress distributions for 
which the 10 resultants and couples have zero values. 

In the second place, quantities may be arbitrarily chosen to 
represent the state of deformation in the plate or shell. In this 
event, it may happen that the 10 resultants and moments defined by 
equations (33) are not sufficient to describe the corresponding 
state of stress, and additional quantities ("higher moments" of the 
stresses) may occur, together with corresponding additional equilibrium, 
equations complementing equations (35) and (38). 

It is evident 'that in various theories of plates and shells, 
developed along such lines, the final formulations may differ in 
consequence of differences in basic ass ump tions . 

It is felt- that neither of the foregoing two points of view has 
heretofore been adequately taken inter account^"" In this report "a 
definite formulation of the theory with reference to the first point - 
of view is not obtained. However, results are arrived at which are 
believed to be definite with reference to the second aspect of the 
general problem. 


CONVENTIONAL ASSUMPTIONS 


It may be expected that a state of small bending and stretching 
of a plate or shell Is described in the first approximation by the 
following formulas for the displacement components occurring in 
equations ( 8 ) and (9): 


(& 1 j £ 

) = u i(^ij £2) + £ u i*( s i* 6 2) 

( £ 

) = u 2( 5 1> §2) + Su2*( 6 1^ *2) 


W (*1j £) = v(t l3 ig) 
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The components of strain are determined in terms of the five displacement 
variables u-j_* u 2 * u 2 * * and w "by introducing equation (4l) into 

equation ( 9 ). Equations (13)* (l 8 )* and ( 19 ) then serve to express the 
components of stress in terms of the displacement variables. If these 
results are introduced into equations ( 33 ) * the resultant 10 equations* 
in addition to the 5 equilibrium equations* comprise 15 equations 
involving the 10 resultants and coupleB and the 5 displacement variables. 

The success of the customary procedure depends on the use of the 
following argument. An order-of-magnitude consideration of the 
equilibrium equations for the stresses shows that* unless the surface 
loads are highly concentrated* the transverse normal stress 0 £ is in 

general of smaller order of magnitude than the stresses CTj_ and cr 2 . 

In consequence of this fact it is conventional to neglect the term 
involving a £ in equations (14) . 

The third of equations (9a) shows that the assumption of equation (4l) 
implies the assumption of = 0. However* a theory which includes the 
two hypotheses cx£ = 0 and = 0 would* in particular, fail to lead to 

correct results in the special case of a flat plate subjected to a state 
of homogeneous bending and stretching* for which problem the exact 
solution is easily established. This difficulty is usually avoided by 
neglecting cr^ in the stress— strain relation of equation (l4a) and by 

then determining €£ from, the resultant equation. To remove the 

resultant inconsistency* it would then be necessary to correct the 
original expression for W by the addition of terms which are linear 
and quadratic in £ . If no boundary layers of width of the order of 
the thickness h are present* these additional terms are found to be 
small in comparison with the leading term w. Thus* to obtain a first 
approximation theory the additional terms may be omitted in intro- 
ducing V into the expressions for the strains 6 -j_* e^* and 7 q 2 * 

If the assumptions of equations (4l) are introduced into the last 
two of equations ( 9 b) there follows 


7 ^ 1 + + a l 



(42) 


It has been customary to assume further that the transverse shear 
stresses t-j_£ and t ^ are also of smaller magnitude than c r 1 * o^* 

and t 2_ 2 . If these stresses are neglected in the stress— strain relations* 
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two of equations ( 13 ) "become 7 ^ = 7 2 £ - 0* With these relatione, 
equations (h2) determine the "bending terms uj_ r and U 2 f as functions 
of U]_, u 2 , and w as follows: ... 


ui‘ 


1 dv U 1 
+ % 


U 2 * 


_ 1 _ dw_ ^2 

a 2 Sl 2 e 2 


( 43 ) 


It may be mentioned that the assumptions e £ = 7 ^ = 7 ^ =0 at 

all points are equivalent to requiring that straight lines which were 
originally normal' to the undeformed middle surface remain straight lines 
after deformation, remain normal to the deformed middle surface, and 
suffer no extension. 


"With these simplifications, the expressions for Qq and Q 2 in 
equations (33) can no longer be retained and there results a set of 
13 equations involving 13 unknown quantities, Uj T and u 2 * being 
eliminated by means of equations (43). 

If it is assumed that cf£ is negligibly small, there are obtained 

from equations ( 13 ), (14), and ( 33 ) the following equations relating the 
stress resultants and couples to the components of strain: 
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Equations (44) are listed as a basis of an outline of certain 
procedures which have been used to obtain relations between the 
stress resultants and couples and the displacement variables u^, ug, 
and w. 


LOVE’S FIRST APPROXIMATION 


.If the ratio R is neglected in comparison with unity in 
equations (44), as well as in equations ( 9 ) which define the strain 
components, the following stress— strain relations are obtained: 


%i - + V6 2°) 

*22 - p^2 ( e 2° * V ‘l°) 

N*j_2 = = Gh7 i 


%1 = + W 2) 

+ VBl ^ 

%2 = %L - O t| t 


> 


(45a) 


(45b) 
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where 


€ 0 „ JL_ ^1 u 2 ^1 w 

1 " a x oqo^ ai 2 


6 0_1^£ , U 1 dog 1 w 
2 02 81 2 0102 81 x e 2 


y 0 _ “2 8 /V\ ai a AlN 

12 «i ^Vo2/ + cc2 ^Vov 


K - 1_ hu l , V_ 52l 
1 a-L a6 x o 1 a 2 ai 2 


(46) 


1 8 u 2 t hi 1 8 ctg 

2 og at 2 “i^ ai x 


= 52 + EL -Z-fin!) 

ct-l aii \02 ) 02 as 2 \ a l / 


and. where uj T and u 2 ’ are defined by equations (43) . 

The system of equations (45), (46), (43), (35), and ( 38 ) was firBt 
given by Love (reference 4, page 531 ; see also reference 5 ) and has 
been used as the basis of many studies of specific problems with 
regard to flat plates, cylindrical shells, and shells of revolution. 

It is generally believed that this formulation of the problem 
contains all the essential facts necessary for the treatment of thin 
shells, as long as special conditions do not require inclusion of the 
effect of transverse shear and normal stresses. 

It is well known that within the framework of this approximate 
theory fewer boundary conditions can be satisfied than is expected. 

It has been shown in earlier papers (references 1 and 2) that for flat 
plates this difficulty is resolved if transverse shear deformation is 
taken into account in an appropriate way. In this report it will be 
shown that the same is true with regard to the analysis of shells. 


MODIFICATION OF LOVE’S FIRST APPROXIMATION 


A number of writers (references 6 , 7 , 8 , and 9 ) have modified 
the foregoing first approximation by not neglecting the ratio £ /R in 
comparison with unity in equations (9) and (44) but still retaining the 
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assumptions of 1 equations (4l) and (43). In particular, if terms 
involving powers of h through the third are retained, and terms 
involving higher powers of h are neglected, inequat ions (44), 
equations (45) are replaced by the following forms (reference 9 ) : 


N 11 = 


Eh 

1 - V r ‘ 


N^ = Gh 


N 2i = Gh 


(<i° * v S2 o) _ - Ei" 

r 12° - 1S(^, ~ ^)(" 12 - ^f-)j 
Ml1 = izf- [( B1 + VB2 ) " & _ jfe) €1 °] 

** ~ aaC?! v 2 ) [(“ s + V,tl ) “ (£ _ ^) 2 °] 
m 12 = G $||t ~ ($£ “ ^) P 1°] 


^21 = G 12 T (r£ Rj) P 2 0 


(47) 


In addition to the quantities defined in equations (46), the 
quantities Pi® and P 2 ® are defined by the equations 



1 chig U 1 3o-i 

Oi ^ ^ 




(48a) 
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and k-j _2 is defined by the equation 


Bo° Bn 0 

2b 12 - T + |2- + U- 


(48b) 


It is noted that the equation 


J3i° + 0 2 ° » ^ 12 ° 


(49) 


1 is satisfied by virtue of equations (48a) and the third of equations (46). 

It is seen that equations (47) differ from equations (45) in. that 
certain terms of order h3 are added to the original expressions. 

However ,’ it may he said that the additional terms which appear in 
equations (47) cannot he expected to he of the same form as terms of the 
same order of magnitude which would he introduced if the simplifying 
assumptions = 7 ^£ = 72 £ = cr^ = 0 were replaced hy more flexible 

assumptions. 

It may he noted that all the additional terms in equations (47) 
disappear when the principal radii of curvature of the middle surface 
are equal, that Is, in the cases of flat plates and spherical shells. 


LOVE’S SECOND APPROXIMATION 


A second approximation, given hy Love (reference 4, page 533 ) > 
introduces three types of corrections to, his first formulation. Love 
states that such modifications are unnecessary unless the flexural 
strains £k-j_, £ 1 * 2 , arbt ^ £ T are large I 11 comparison with the extensions! 
strains eq°, € 9 °^ and 7 qp°. 

As a first modification, the transverse displacement W Is 
expressed in a more flexible form. 


0 = v (*l* * 2 ) + ^( l l* hj 0 
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according to which the direct strains, given by equations ( 9 a), 
■become 


” r +V^i ^ 1 - 0 + Clt:L + |[) 


'2 


“ i tVa s ( 62 ° + E " 2 + 


(h9&) 


- _ <& 
e Z " w 


The second modification consists in not completely neglecting the 

ratio £ /R -with respect to unity hut writing 1 ■ — Z 1 — Z /H-, . 

1 + £/Ri * 

However, in expanding the first two of equations (49a) Love assumes 

t ft Z V 

that the quantities ^ 6 and ^ ^ can he neglected. In this way, 
equations (49a) are approximated in the form 



(49h) 


The term v is considered as a small correction term. To obtain a 
first approximation to its value, the transverse normal stress cr£ 
is neglected in equation ( 19 ) and gj + *2 ^©Placed hy its first 

approximation + C + ** 2 ) • '^ 1 © resu l'* : then becomes 





V£E 

(l - v ) e ^ 


( 6 1 ° + e 2 °) + £( k 1 4 * 2 ) 


(50) 
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The correction w is then determined by integration, with the 
convention v(i*i, °) = 0 l 11 tha form 


w 


VfE 


(1 - v)E^ 


k< 


e i° + 


6 2 °) + i £ 2 ( k 1 + k 2 ) 


( 51 ) 


If this result is introduced into the first two of equations (191) 


a-nfl quantities of the form 



are again neglected, there follows 


*1 = g i° + (s - §f) K i - i (i - v)^- ^ 2 ( k i + “ 2 ) 

0 + U - - 1 ^2 £ 2 (K n + K ) 

V 2 2 (1 -v)E£ ^ 1 + 


e 2 “ e 2 


( 52 ) 


The third modification consists in not neglecting in 

equations (14) . In order to obtain a first approximation to the value 
of cr£, the third of the equilibrium equations (ll) is used. It is 

supposed that no body force is present and that the transverse 

shear stresses are negligible. Also, to this approximation the 
ratios £/R can be neglected, and the equation becomes 


^ ff £ _ SJk + 

ar ~t + 


£2 

e 2 


If cr^ and cr 2 are replaced by their first approximations and the 
extensional strains are neglected, then there follows 


I? = I^[4[ (*1 + VK s) + 4( K 2 + VK l) 
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In the case of vanishing normal surface loads, 5 2 , — 

Love obtains by integration 


_ _ 1 E ft 2 _ g\ Aj + VKg 
2 ! _ V 2U 5 /\ E 2 


+ 


K 2 


+ VK l\ 

e 2 ; 


(53) 


If the expressions for C]_ and. cr 2 obtained, by introducing 

equations ( 52 ) and ( 53 ) into equations (l4) are in turn introduced 
into the first two of equations ( 33 ) and if again quantities of the 

form ^ are neglected , there follows finally 



together with an analogous expression for N 22 . 


A comparison of-equations (54) with the first of equations (47) shows 
that equation (54) includes the terms present in the former solution 


except for the term 



which was neglected in the present analysis. 


More important j however , Is the fact that new terms, of the same order as 
the correction terms present in equation (47) , are introduced in 
equation ( 54) in consequence of the partial inclusion of the effect of 
the transverse normal stress d£. Furthermore, these terms do not 

vanish when E-j_ = E 2 . Still no account has been taken of the possible 
effect of the transverse shear stresses t <, and r ot , . 


The expressions obtained, in the present procedure, for the 
remaining quantities listed in equations (47), are equivalent to the 

e 0 *y0 qO 

results of neglecting terms of the form and ■%- in comparison 

with the flexural terms in equations (47) . That is, no additional 
terms due to the effect of cr^ are introduced Into the remaining 

quantities by this procedure. 
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It may be mentioned that the results of the developments of this 
section reduce to those presented In Love*s treatise in the isotropic 
case = V, E^ = E, which Is the only case treated there. 

BASSET’S THEORY 


In this section there is outlined an analysis given hy Basset 
(reference 10) which, in the opinion of the present authors, has not 
received as much attention as It deserves. The reason for this may 
perhaps he found in the fact that Basset’s wort is difficult to read 
and that the notation employed is somewhat complicated and, from 
modern standards , somewhat unsystematic. 

Basset begins his analysis with the assumption that the stress, 
strain, and displacement components in a shell can he expanded in 
series of powers of £. Thus there may he written, for example. 


u i(^ij ^2-> £) = u i(£i* £2) + £ui* (§1, 62) + \ f; 2 *! 11 (iij 12)+ • • • ( 55 ) 

where 

ui(&i, 12) = *i(si, £2, 0) 

U1 ’( ei ' 52 ) = i i ] c= o 


and so forth. 

The derivation is based on the use of the principle of mi nimum 
potential energy (equations (25) to (28)). The strain-energy function P, 
as given hy equation (27), Is expanded in powers of £ so that the first 
terms are of the form 

P = 2 |i I v2 ' 1 ( 6l + ^ 6l * + 2 + * • *) 2 + ( 6 2 + ^2’ 

( 57 ) 
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where ei is written for e^lj, ^ 2 > °) > e l*> for — € — — - 2 — - ; and 

bo forth. Only such terms are retained aa will lead to tarns of 
order h 3 or lower after the £ —integration of equation (25b) is 
carried out, terms involving the basic quantities u-^, u 2 , and w being 

considered as of zero order in h. 

To express the primed quantities in terms of unprimed quantities, 
use is made of the strainr-displaeement relations of equations (9b), the 
last two of which can be written in the form 


Cftbl •] 'vrj U-. 

bT = ai(l + "S/R a J) dIJ + ST^T 


du P 

sr = r 2£ - 


aw + u 2 




a^ ^ ag^l + £/Rgj ^2 Eg + c 

Thus, setting £ = 0 , there follows 

u . t = 7 , t _ _1_ aw_ JJl 

ug T = 7gt ^ 

^ ^ ag d £ 2 E 2 


( 58 ) 


( 59 ) 


In this way the quantities uq* and u 2 ’ are expressed in terms of the 
displacement functions U]_, ug, and w on the middle surface and in 
terms of the transverse shearing strains on the middle surface. The 
third of equations (9a) gives the result 


w’ =* (60a) 

This result can be expressed in a different form if equation (19) is 
used to express in terms of € g , and cf£ as follows: 


w f 




(60b) 


By further £ -differentiation, all other primed quantities can be 
expressed in terms of- the values of u^. Ug, w, and fa +- e 2 ) on the 
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middle surface; such expressions also Involve the values of 7-^, 7 2 £ » 
ftriri cr^ and their £ —derivatives on the middle surface. The remaining 

three stress— strain relations of equations (9) then permit the 
determination of the coefficients in the power-series expansions of 
the strains ^2* arj & 7 io in terms of the j 3 ame quantities. 

It is then assumed that in the expansion 

^(^Ij hj 5) = h) + S^'(ii, 62) + I £ 2ff g” (*i> *2) + • • • 


the leading term is at least of order h^ and the coefficient of £ 
is at least of order h. Next it is pointed out that the leading 
terms in the expressions for the couples Mjj_, Mpp, M, p , and 

are at least of order h. 3 . Qn the "basis of this fact it 1 b stated 
that, if the external moments m^ and ioq are absent, the equilibrium 
equations (38) imply that the transverse shear re stilt ants Qq_ and Q2 

must be at least of order h .3 and hence that the terms of lowest 
order in the transverse shearing strains 7 j_£ and 7 2 £ mast be at 
least quadratic functions of h and X since such functions when 
integrated over a section give rise to quantities of order h 3 . 


The fact may be pointed out here that this argument may break: 
down in those cases where an appreciable change in the magnitude of 
a couple may occur. over a distance of the order of magnitude of the 
thickness h of the plate or shell or over a distance of order \/ah 
where a is a representative dimension of the shell. For if such a 
change takes place over a distance l in the i^-direction it follows 


that 4 ^- is of the order ~ M in the region considered. It is well 

known that for plates l may be of order h, while for cylindrical and 
spherical shells the distance may be of order \|ah, where a is the 
radius of the circles of curvature, in problems of usual occurrence. 


If such cases are excluded, it is found that the contributions 
of and to the strain-energy function of equation (57) give 

rise to terms of order bP or higher and hence the transverse shear 
effect may be neglected if only terms of order h 3 are to be retained. 

"With the assumptions noted, the variational equation (28) leads to 
expressions for the stress resultants and couples which include al 1 
third— order corrections to Love*s first approximation, equations (15), 
which are consistent with the assumed orden-of-magnitude relations 
involving 7^3 7 2 £, and cr£. In particular, the corrections introduced 
"by .equations (lj-7) ere included, as well as the additional corrections 
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introduced by Love’s second approximation, and the remaining third- 
order corrections which are absent in the latter approximation are 
introduced . 

Basset’s analysis was carried out only for isotropic circular 
cylindrical and spherical ^shells "but could equally well he adapted to 
the more general theory. 


TOWT 7 ,* THEORY 


An alternate system of equations for small deflections of isotropic 
shells has been derived by E, Trefftz (reference ll) by means of the 
minimum principle for the stresses. Trefftz begins by writing 


a l = + S^l 1 


0g = +■ (JCg® ^ 


t 12 " T 1 2° + ? t 12' 


(61) 


where the functions 0® and cr* are expressed in terms of resultants 
and couples by means of equations ( 33 ) • Equations ( 6 l) are introduced 
into expressions (30) and (32) for the complementary energy, in which 
all terms containing cr^, and are omitted. 

The resultant-expression for the complementary energy expressed in 
terms of the N’s and M’s is minimized subject to the restrictions 
imposed by the five equilibrium equations (35) and (38). The equilibrium 
equations are taicen into account by means of the Lagrangian multiplier 
method, the multipliers being identified with appropriate displacement- 
components . 

When termB of order £/R are neglected, Trefftz’ results agree 
with Love’s first approximation as given by equations (k$) and ( 46 ). 

When terms with £/R are retained a new system of equations is obtained. 

With regard to this set of equations, however, the following point 
may be made. For certain exact solutions for circular rings it is found 
that, as h/R increases, the linear displacement expressions corre- 
sponding to equations ( 4 l) are much more nearly correct than the linear 
stress expressions corresponding to equations ( 6 l). It is noted, 
however, that this distinction disappears for the special case of the 
flat plate. 
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The effect of transverse shear and normal stress In shell theory 
might he taken into account hy calculating hy means of the three- 
dimensional equations of equilibrium (ll) the values of the stresses cr«, 9 

t and t 0 -which, correspond to the values of or, , cr_, and as ^ 

- L b _ 1 / 2 * 12 

given hy equations (6l) . 

It appears that this procedure would lead to rather complicated 
expressions for cr^, and t 2 £ when R 1 and K 2 are finite. A 

class of non h omogene ous shells for which this complication can he avoided 
has been studied in reference 12. 

In earlier papers (references 1 and 2) one of the present 
authors has treated the problem of the flat plate in this manner and 
has thereby established the importance of the role which is played by 
the deformations caused by the stresses T l£» r 2 t cer 'fc8-i n 

problems of plate theory. It may be remarked that the use of the 
particular minimum principle employed in that reference is not essential 
to the analysis of the effects of transverse stresses. Similar results 
can be obtained in a number of ways, for example, as is shown in a later 
section of this report, by use of the mini muni principle for the 
displacements . 

It may be added that the system of equations obtained by the 
application of the minimum principle for the stresses in the manner 
Just outlined, while furnishing approximations for the stresses 
themselves, determines certain weighted averages (taken over the 
thickness) of the displacements, rather than the displacements themselves. 
This is shown in reference 2 with regard to flat plates. It will be 
shown in the present report that Just the reverse is true when the 
minimum principle for the displacements is used to obtain a system of 
two-dimensional equations for the theory of shells. 


THEORY OE SYNGE AND CHIEN 


A theory of the finite deflections of isotropic shells has been 
developed in a series of recent papers by J. L. Synge (reference 13) 
and 'W. C. Chien (references 13 and l 4 ) . The basis of this theory, 
which Is of great generality, appears in the case of small deflections 
to reduce to the following considerations. 

All stresses and displacements are expressed as power series 
in £ and the equilibrium equations and stress— strain relations are 
written accordingly as equalities between power series in £ . By 
equating the coefficients of respective powers of £ in each of these 
equations, an infinite system of simultaneous differential equations 
in the independent variables £-j_ and t 2 is obtained. 
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Next a representative length a on the middle surface of the shell 
is introduced and each coefficient of each power series in £ is 
considered to he developed in powers of the ratio = h/a. Various 
choices of the lowest exponent (which may he negative) of- q are made, 
such choices leading to different formulations of shell theory. 

Finally, only the leading terms of the respective developments are 
retained . 

The fact must be stressed that this interpretation of the Synge— 
Chien procedure is tentative and that a more detailed explanation of 
this work in terms of the notions of the present report is believed to 
be desirable. It may be remarked that the formulation of boundary 
conditions appropriate to the various theories evolved by this procedure 
does not appear to be incorporated in the work of Synge and Chien. 


^CORPORATION OF TRANSVERSE SHEAR AND NORMAL STRESS EFFECTS 
WITHOUT USE OF A MINIMUM ERINCIFLE 


In this section there is outlined a possible extension of the 
classical general methods which takes into account a first approximation 
to the effects of transverse shear and normal stresses. 

Again there are assumed for the displacements the expressions 


U l( 5 l* 0 = Ui(6i, 1 2 ) + £ u l’(*l> 6 2 ) 
u 2(*1j ^ 2 > 0 = *2) + £ u 2*( S 1j £2) 

w (*1j * 2 * 0 “ w (*l* h) 


> 


( 62 ) 


It is hot, however, assumed that the transverse shear strains are 
negligible. Thus equations (^ 3 ) no longer hold. In their place 
now there are the definitions of Qj_ and Q2 as given in equations (33). 
To a first approximation these equations take the form 
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Qi 


Qe 


p/2 

<% | it 

J-h/2 


1 


Ph/S 




-h/2 




(63) 


where the factors 1 + £/R are replaced "by unity. With the same 


appr oxj mat ions 

, expressions 

for 


and 

equations (9h) 

in the form 





7 lt 

1 

8w 

U 1 


ai 


"SI 



= X 

og 

8w 

5I 2 

u 2 

Thus } in place 

of equations 

(33) 

now 

there 


ui' 

+ k 

cftr 

_2L , 
E 1 



1 

8w 

u 2 


u 2 

“2 

d| 2 

e 2 ' 


(64) 


( 65 ) 


The next step consists in introducing the values of e 2j 

and t-j _2 which correspond to equations (4l) into the firBt of 

equations ( 13 ) and into equations (l4)j to determine corresponding 
expressions for cr-^ cr 2j and 7 -j P and in then substituting the results 

into equations (33)., to determine corresponding expressions for the 
stress resultants and couples . In doing this., again the factors 1 + £ /R 
are replaced by unity in calculating terms which are in addition to 
those which appear in equations (47) . If the right-hand members of 

equations (47) are denoted by W pp® 3 and so forth and the parameter 
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VfE 

j* = k_ 


v* = 


(1 - V)E^ 


( 66 ) 


which appears in equations (15) and (16) is used, the resultant 
equations take the form 


lh/2 

N 11 “ %° + V* ' <t £ 45 

L — h/2 


Ngg = W 22 ° + v * 


»12 - V 3 


%L - %L° 


nh/2 


U- h/2 


d£ &£ 


(67a) 


Hh/2 

M u = M 11° + V* I cr^ d£ 

U~h/2 

ph/2 

H> 2 “ ^22° + V* / d £ S d£ 

U— h/2 

^ = M 12 ° 

^21 = %L° 


(67b) 
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Equations (67) and (63) , together with the . 5 equilibrium equations ( 35 ) 
and (38), comprise 15 equations,, involving a total of 17 unknown quantities, 
namely: 10 stress resultants and couples, the 2 integrals involving cr^ 

in equations (67), and the 5 displacement variables of equations (62). 

Two additional equations may he obtained by requiring that the 
results of integrating equation (19), and the product of £ and the 
two members of equation (19), be satisfied. Since the third assumption 
of equations (62) implies the vanishing of the transverse direct 
strain 6£, these equations take the form 



Equations (67), (68), and (63) serve to determine the conventional 
10 stress resultants and couples, as well as the 2 auxiliary resultants, 
in terms of the 5 displacement functions. The introduction of these 
results into the five equilibrium equations (35) and (38) then leads 
to a set of five differential equations in the five displacement 
functions. The solution of the given set of equations is then apparently 
determinate if these five functions are prescribed along the boundary of 
the plate or shell. Alternatively it may be expected that along an 
edge £q = Constant the five quantities N.^, N 12 , , M, p , and 

may be independently prescribed, as statical considerations would 
require. If this is indeed the case the classical difficulties first 
investigated by Kirchoff would be resolved. 


A new difficulty may, however, be noted. 3 h the special case of 
pure bending of a flat plate, by couples distributed uniformly along 
the boundary, the known exact solution is derived from equations (67b) 
only if the corrective terms involving cr £ are absent, that is, 

if Mjj = M-, i° and M^ 2 = ^ 22 ° • ^ equations (68) are to be 


satisfied it is readily verified that the corrective terms will be 
absent only in the special case when the prescribed uniform edge 
couples satisfy the relation Mj_q + Mq2 = 0 . The presence of this 
difficulty may be explained as follows. The assumption W = w^£q, £ 2 ^ 

implies the assumption €£ = 0, and hence (if the relevant stress- 
strain relation is retained) there follows also 
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Thus a transverse normal stress effect Is automatically introduced, 

in this formulation. In connection with effects due to the 

stress cr^ + o 2 . This coupling is avoidedin the classical theory hy 

disregarding the third stress— strain relation and taking cr^ = 0 

identically, even though this in general contradicts the 
assumption = O which is also part of the classical theory. A 

similar procedure might "be adopted here, wherein equations (68) would 
he disregarded and an expression for such as that used hy Love 

(equation ( 53 )) in his second approximation would he used. The validity 
of such a procedure would, however, he open to question. 


It is clear that the difficulty Is not present in the special 
ortho tropic case when is negligibly small. 


A second possible procedure consists in replacing the 

dW 

assumption = 0 hy a more general type of assumption under which 


and c are not so restricted. In this connection, however, it may 

he remarked thatrd.f the assumptions of equations (6 2) were replaced 
hy more flexible assumptions, in which more than five displacement 
functions were involved, the present procedure would not supply in a 
rational way the additional relations needed for the determination of 
the additional functions. 


n 


To clarify the entire situation and, in particular, to investigate 
in a rational way the manner in which boundary conditions may he ' 
imposed, the formulation of the problem is considered next from the 
point of view of the principle of minimum potential energy. 


APPLICATION OF TEE PRINCIPLE OF MINIMUM POTENTIAL ENERGY 


In place of the assumptions of equations (62), here the following 
more flexible forms are taken: 

U l(^lJ 0 = Ig) + Suq'^q, £g) + \ *g) 

u 2(^1j 0 * u 2_(^1j ^ 2 ) + ^ u 2 f (^lJ ^ 2 ) + 2 ^ 2u 2 M ( S l* ^ 2 ) 

w (^i> h* 0 -v(*i, h) + &'(*!* h) +|^ 2v?, ( l i. * 2 ) 


r (69) 
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and a system of equations determining the nine displacement functions 
present is derived "by an application of the principle of TniniTmjm 
potential energy (equations (25) to (28) ) . 

A shell is considered "whose boundary edges lie in coordinate 
surfaces 1 ^. = Constant and £2 = Constant, and there is first 

calculated the potential energy associated with the load Bystem. 

According to equations (25a), 


«7 = — 


(? 1 + U 1 + + P 2 + U 2 + + ^ W +)( 1 + ^)( 1 + 2 ^ 2 ) 

( P 1- TJ 1— + P 2- IJ 2- + ” 2R^)j a l a 2 d ^l 


- 


nn/a 


llJ-h/2 


+ TafcUt + T n? w) (l + Jj.) 4£ 


Ofc d ^t 


= rt, + ir 


Is ^ *le 


( 70 ) 


The first (double) 

integral ti-^ b 

represents the energy associated with 

the surface loads 

p l+* P 2+* and 

* + 

acting on the 

surface 

t = «* and 
* 2 

the surface loads 
The abbreviations 

V ± _, P 2 _i and 

q_l_ 

acting on the 

surface 

t - _ fc 
l 2* 


^1+ = U l(«H «£> + |) U l- - U l(h, - I) 


and so forth have been introduced. The second integral in 

equation (70) represents the energy associated with the edge stresses. 
Here the subscripts n and t refer to “the normal and tangential 
directions on the boundary faces; thus, on a face £^_ = Constant, 
n = 1, t = 2. 

If the assumptions of equations (69) are intrQduced into the 
expression for jt-^ 0 the coefficients of the displacement functions 
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are found to involve quantities which may be denoted as follovs: 


p l 

= P l + < 

^ + ^[)( 1 

+ A) 

+ p i- 

C 1 ^iX 1 

. h_ 

2Rr 

■) 

P 2 

= P 2+ 

v 1 + AX 

+ A) 

+ P2_ 

( x ■ ^X 1 

_ h 

2Er 

•) ' 

q 

■ < l+( 1 

*-AX^ 

A) + 

q-(l- 

41 X 1 2 * 2 ) 


“i = “ 1 


+ AX 1 + 

A)~ 

Pl-(l 

I 

Hsf 

1 

A)_ 


*2 = i[ 

P 2+ (- 1 

+ wX 1 + 

A) _ 

P2-( 1 

AX 1 

A) 

> 

n = 1^ 

( X + AX + A, 

)-*-( 

1 k. 

L 2R- 

lX 1 - A) 




( 71 a) 


(TO) 


With these abbreviations, there follows 


Zs 


- +-52^2 + qw +- m-jUi* + + nw* 

+ 1 h 2 pi^i M + | b 2 P 2 V + £ h 2 qv”^cc 1 a 2 d^ d £ 2 ( 72 ) 


where the integration is taken over the area of the middle surface of 
the plate or shell. 

If the assumptions of equation (69) are introduced into the 
expression for and if the £— integration is carried out, the 

coefficients of the displacement functions involve quantities which may 
be denoted as follows: 




Of the 18 "resultants" so defined, corresponding to edges - Constant and £ 2 = Constant, it is 
seen that 10 are the familiar stress resultants and couples of equations (33) j while the remaining 
8 quantities are new. 

With the notation of equations (73 ) , there follows 
He m " J (^rna u n + *WV + *mPn" + ^nt^t + *ktV + ^ntut' 1 + Qd.w + SqV 1 + T n w» , )a t d£ t (7*0 


where the line integral is taken along the coordinate curves which form the complete "boundary of the 
middle surface. The fact that edge values of the resultants are involved in equations (7*0 is 
indicated "by the use of "bars. ' 


H 
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The internal energy Strain iB gi Yen by the volume integral 


nn 


lh/2 


^strain ~ 


\JU 


[d-h/2 


+ -4) aS 


dg 2 dg 2 


( 75 ) 


■where P represents the strain energy per uniir volume and, is defined by 
equation (27). Thus, the variational principle of equation ( 28 ) is of the 
form 


Blt strain + Q *l B + b *le = 0 


( 76 ) 


■where 


^strain - 


)h/2 


J-h/2 


(sp) ai a 2 (l + ^)( 1 + 


dg-^_ d£ 2 


( 77 ) 


If-use is made of equations ( 2 6), the expression for SP can be 
put in the form 


5P 


dP o .. SP 

8^! 61 de 2 


&€2 "f — • 





= cr 1 Se 1 + a 2 Se2 +-cr^8e^ + T 12 5 ^12 +" T l£ 5 ^l£ +: T 2£ s 7 2 f; (78) 


where the stress and strain components are to be considered as functions 
of the nine displacement functions appearing in equations (69). 

First there is calculated the contribution of the first term 
to the variation . If the assumptions of— equations (69) are 

introduced into the first of equations (9a), there is obtained 
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^3 


e l 


a i(l + £/%) 


(Ifi + E sir + 2 £2 %f) + 4 + 5us ' + ^ ?Su2 ") 


+ ^(w + ^w f + i ^ 2 v'*) 


Thus, the result of replacing BP h y cr^Be^ in the right-hand member of 
equation (77) is of the form 


th/2 


l-h/2 


•C 1 + 4) 


SBut BBut * i „p SBut ** 

^SF 1 + « + 1 ? -si?- 


86 - 


^^(Sug + ^Bup* + 7? £ 2 Bu2 ,T ) + ^“(&v + £bw’ 


+ 




dddi x dl 2 


If the £— integration is carried out and the notation of equations ( 73 ) is 
introduced, this expression takes the form 


N- 


11 


(- S? * ft * ft -■ *¥ -) 


88ui Tt Sag 




5^2 8v „ 


46 ! 4§ 2 


where the integration is taken over the middle surface of the plate or , 
shell. 

The contributions of the remaining terms in equation (78) to the 
value of 6jr B train. can calculated in a similar way. In calculating 

the contribution of the term t -j p Sy^g ii* is convenient to make use of 
the identities of equations (7)* 
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The contribution of the term cr^Se^ involves two new "resultants, " 

not defined by equations (73), for which the following notations are 
adopted : 



( 79 ) 


In the resultant expression for ^gtrain the terms involving 
derivatives of the independent variations may be integrated by parts, 
this procedure leading to an equivalent expression for train 

consisting of the sum of- a line integral taken around the boundary of the 
middle surface and a new double integral involving the independent 
variations in a linear way. 


Details of the calculation are omitted, and only the form is 
presented that - is assumed by equation (76) when the separate integrals 
are combined and the coefficients of the independent-variations are 
collected: 
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dap N-| p Sa]_^22 

_ ^ t ± 3t 2 


+ | + W 2 i - N 1X |2L + Q2 + ai 02P2 

b l 


dai 
dig 


Bug 


&v 


^ + ^_(^ + | £ ^ a2 + aiaaa - 


+ Sogu. + Sa^i + ^ |g _ ^ |5£_ Qiaiaa _ aia2mi 


da^M^g ^ a l'^22 


dig 

dag 


da' 


an + ~htr + ^ sn " * a ata " Q£ ° ia2 + <Wa 


5u^* 


Bu^ 


* -'7^ - (ffj 1 + - "“l-s * “l“S D ] 6 ''’ 

Sf 11 * ^ - ** is - ("1 * %)•« 


+ | a-LOghSi 


Bui ** + 


do2?i2 . ^l p 22 . “n . dccg „ da-, 

"s TT ~ %T 21 ail " 11 Sei 


- ( s 2 + + 1 “1 “2^12] ^a" + pgf ^ 1 + 

- (gj 1 + ^ £ ) a l a 2 “ &WS + | ttia^ll&w" 


dlj. dl 2 


-m — ^nn)®' u n + (®nt ~ -^nt;) ^ u t + (^1 ~ On)^ 

+ (*xm- *hn)&V + (%.t ~ %t)&*t r + (^ - S^Sv* 

+ (^nn ~ p nn)^ u n ,f + (-^nt “ + (^n ~ ^n}^ w ’ , ^|‘ a t 


= 0 


(60) 
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Since the -variations appearing in this equation are entirely 
arhitrary in the interior of the middle surface, the coefficients of 
the nine variations in the double integral must each vanish identically. 
In this way nine differential equations are obtained, the first five of 
which are the conventional equilibrium equations (35) and- (38), except 
for the interpretation of the loading terms, and the last four of which 
are new. These additional equations correspond to the introduction of 
the four new displacement functions w 1 , w M , u^**, and u^**. 

The corresponding nine boundary conditions to be prescribed along 
the edges of the plate or shell are obtained as a consequence of - the 
independent vanishing of the nine terms in the line integral of 
equations (80) . Thus, along a boundary = Constant, the first term 
of the line integral becomes — N q q ) 5 u j . The vanishing of this term 

is assured if the displacement = Ui(iq, ip, 0) is prescribed along 
this boundary, since then Buq = 0 . Alternatively, the term will vanish 
if the resultant- Nq q is required to take on a prescribed value N-| q 
along this boundary. Considering the other terms in the same way, it is 
seen that - the following nine boundary conditions along a boundary 5 q = C 
are consistent with the displacement assumpt ions of equations (69): 


Nu = W 11 

or 

= 03.(0, 62, 0 ) 

prescribed 

%2 = %2 

or 

u 2 - U2(C, i 2 , 0) 

prescribed 

11 

or 

W - W(C, l s , 0 ) 

prescribed 

& 

H 

11 

f 

or 

au^c, tg, 0) 
1 3? 

prescribed 

II 

CM 

s 1 

or 

au 2 (c, ig, 0 ) 
2 35 

prescribed 

Si = s x ■ 

or 

... awp, 0) 
dS 

prescribed 


Tl= Tq_ 


Pn = *11 


or 


or 


or 


w 


it _ 


tPw(C, j 2 , 0 ) 

a ? 2 


Ul « a^i(c, h, Q) 
1 sp 

„ a% 2 (c, i 2 , 0) 

w — 


prescribed 


prescribed 


'(8lb) 


P 12 = *12 


prescribed 
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In particular, it is seen that if w 1 = w ,T = u^” = vlq ' 1 = 0 

everywhere, and hence equations (69) are specialized to the classical 
assumptions of equations ( 4 l) , the five "boundary conditions which remain 
are those of equations (8la) . Thus, for example, at a free edge of a 
plate the twisting couple M-| P and the shearing force Qj can "be 


required to vanish independently. These results are in contrast with 
the corresponding situation in the classical theory, where Kirchoff 
first showed that for a flat plate the physically desirable 
conditions Mj_ 2 = Qi = 0 at a free edge cannob both be satisfied but 


that they must be replaced by the single condition 


Qq + 


8^2 


= 0. 


Basset, Lamb, and others have since shown that in the classical theory 
the resultants Mgp , and Qj_ cannot in general be separately 

prescribed along an edge of a shell. (See reference 4 , pp. 536 and 537 .) 


The next step in the analysis consists in expressing the resultants 
of equations (73) (79) in terms of the nine displacement functions. 

The expressions for the relevant stress components are obtained by first 
introducting equations (69) into equations (9) to determine the strain 
components and then substituting those results into the stress— strain 
relations of equations (13) and ( 18 ) . The remainder of the analysis 
consists basically In introducing the expressions for the resultants 
into the equations obtained by setting the nine bracketed expressions 
in equation (80) equal to zero and so obtaining nine differential 
equations in the nine displacement functions. In practice, however, 
the latter part of the analysis may frequently be carried out more 
conveniently In terms of some other appropriate set of nine independent 
quantities. 


Explicit expressions for the "auxiliary resultants" A and B 
of equations ( 79 ) t which do not enter into the boundary conditions, may 
be listed as follows: 


A 


B 


Nn + N ' 


22 



R 2 ) 





2P22\ 
R 2 ) 



In obtaining these expressions, use was made of the third of 
equations ( 12 ) to express cr^ In terms of (cxi + cxo) and <=£. 
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By introducing equations (69) into equations (9) expressions are 
obtained for the strain components which may be written in the form 


61 = r^D^( 5i ° + £s i’ + 1 £2e i") 

62 m Z 7 tK ^° + tes ' * ® £262 ") 

€£ = v’ + £w ,f 

712 " + tfV + 2iapi ") + TTijrM + CP2 ’ + i ^ 2 ”) 

7 ^ = + + 1 'V') 

~ rrtT^ 0 + Cp 2' + 1 ^") 

where the relevant functions are defined by the equations 


6 . 0 _ _1_ , _v 

e 7 f = JL ^1 , w* 

1 a l a l<*2 b^ ^1 


1 


(83) 


tt = i ^1" , V 1 ^1 , V 11 

^ a l S «-i“2 85 2 E 1 


( 84 a) 
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6,0 = 

^ *1 dlq CXjOg 


Ban 


* = JL ^g 1 _ u l* 5a] 


Pi’ - 


a l 51 -l c^ag dig 


(84b) 


Hi 


Hi' = 


_ 1 

Su 2 m 

a l 

56i 

= -1- 

dw _ 

«1 

^ll 

_ _1_ 

dw* , 

ai 

5 Si 


+ Hi 1 


M 


(84c) 


n , ,t = JL &lIL + Iil!l 

41 ai dli + Ei 


and by the additional equations obtained by permuting the subscripts 1 
and 2 in equations (84). 

The calculation of the resultants involves integrals for wl^ich 
there are introduced the following notations 
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together with the analogous expressions obtained by permutation of 
subscripts. The par am eters introduced in- equations (85) possess the 
following expansions in powers of h: 

ai ° = 1 + 5(21*1) + 7(261) 4 ~l(p6i) + * * ' 

-i- - i + f^) 2 + l{A) k + 4(%) 6 - • ■ 

a i" = 1 + K^fe) 2 + 4(261/ + 4(261/ + * • • 
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■With the notation of equations (84) and (86), the expressions for 
the resultants defined "by equations (73) can "be written as follows: 
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The expressions for the nine remaining resultants are. obtained by 
permutation of subscripts. 

Ihe results listed in equations (87a) reduce to those of equations (V?) 
if w* = u^'* = ug” = v 1 * = Oj if up* and U2 r are determined from 
equations (^3), if = 0^ and then if terms involving, powers of h 
greater than three are neglected. 

In this connection it is useful to write the expressions 
for_ Mq, and Ppp in the following Alternative forms: 
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where 


V* = 


E* = 


VgE 

Cl - V)E£ 

(1 - y 2 )Hg 

1 - 2V^V* 


(89) 


The coefficients of v* in equations ( 88 ) represent the 
contributions of - the transverse normal stress cr^ and. are zero 
when cr^ vanishes. In accordance with the remarks made in the 
preceding section, it is evident that if the effects of the transverse 
normal stress are to "be taken into account in a reasonably accurate 
way, the terms w’ and w ,f must not be excluded, in general, unless 
the orthotrophy is such that V£ is negligibly small. 


It appears, however, that the contributions of the terms U ]_* 1 
and U 2 m are in general of small importance. In fact, it is seen 
that these terms appear, principally in combinations of the 

form uq° + u^*' and u 2 ° + ^ u 2 n , such combinations representing 

the average values of the assumed displacements Uq and U 2 over the 
thickness. Thus, it appears that-neglect of the terms Uq ft and u 2 M 
can be largely compensated by interpreting the terms Uq° and u 2 0 

as average values rather than as values assumed at the middle surface. 
The equations herein listed afford the basis of a more detailed 
investigation of the importance of the terms Uq* * and u 2 ,f . 


If these terms are neglected there remain seven displacement 
functions and, consequently, seven boundary conditions are to be 
satisfied along an edge of the plate or shell. Thus two conditions 
in addition to the usual five conditions suggested by statics are to 
be imposed. In particular, if the boundary conditions Involve the 
stresses, the two additional conditions prescribe the integrated first 
and second moments S and T of the edge transverse shearing stresses, 
in addition to the transverse shear resultant Q. If the terms Uq M 
and u 2 M are also retained, the second moments of the edge normal and 
shear stresses in the direction of the middle surface are also to be 
prescribed. 

It is of some interest to investigate in what way the number of 
boundary conditions is reduced when the effect of transverse shear is 
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neglected "by talcing G-£ = w. Eeference to equations (87) tlien shows that 

in order that the resultants Q, S, and T remain finite there must he 
taken 


hi 0 = Hi* = Mi*' = 0 


( 90 ) 


and hence. 


in accordance with 


equations ( 81 c) , there must follow 


u-t* = ^ i- ^r- u 2’ = 

1 ®1 S&j_ ^ ®2 S |2 


U]/* = 


*r sir 


up" = - J- 
^ “2 Sip 


> 


(91a) 


and also 


1 Sw* __ Sw" _ _1_ Sw* _ Sw" _ a 

^l S|q_ S|]_ ^2 Sip Sip 


( 9 lt>) 


The nature of the relevant boundary conditions is determined by 
considering the form assumed by the variation of the expression 7 tj e 

which represents the energy associated with the edge loads. If use 
is made of equations ( 7 l) and (91a), this variation becomes 


5 n e = - ? 


- £ ffc) - £ w §£ 


+ N nt Su t + “a t f^~ 


c 1 SSw N 


+ Qq&w + SjjSw* + T n 5 v" 


Ofc &St 
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The terms involving tangential differentiation can "be integrated by 
parts to give the result _ * . 



Ey considering the results of the application of the principle of 
minimum potential energy, it is concluded that either the quantities 
whose variations appear in equation (92) must be prescribed or the 
quantities which appear as their coefficients must be prescribed edge 
values of relevant resultants. Thus, in particular, if transverse 
shear effects are neglected, the three resultants Q n , N n t, and 

cannot in general be separately prescribed at a boundary, but only 
"effective" shear resultants of the form 


( 93 ) 


N n t - *nt + 

"t 


% = a + xiM 

^ ^ + at 


can be specified. This result is in accordance with known facts first - 
established by Kirchoff for- the flat plate and by Basset (reference 10 ) 
for the circular cylindrical and spherical shell. 


CONCLUDING REMARKS 


From a survey of various systems of equations given in the literature 
for the analysis of small deflections of thin elastic shells, it appears 
that the question concerning the best form of the basic system of 
equations of shell theory has not yet been decided, even in the small- 
deflection theory. 
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The present authors believe that their approach by way of the 
principle of minimum potential energy represents definite progress 
for the following reasons: (l) It is now possible to have a succession 

of two-dimensional theories of varying degrees of exactness, depending 
on the number of terms which are retained in the series of the 
displacement components and (2) the equations which are obtained include 
the possibility of analyzing boundary— layer effects not only when the 
boundary layer is of order \ISh but also when the boundary iayer is of 
order h, where E Is the radius of curvature and h is the thickness. 


Massachusetts Institute of Technology 
Cambridge, Mass., June 23, I9V7 
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